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FOREWORD 

This Indian Standard was adopted by the Bureau of Indian Standards, after the draft finalized by the Hydrometry 
Sectional Committee had been approved by the Water Resources Division Council. 

Unsteady flow models are appropriate for computing stream flow records at locations where: 

a) a single-valued stage-discharge relation does not exist, 

b) backwater affects the discharge under selected or all conditions, 

c) flows are affected by tides, or 

d) it is not possible to gauge the flow using velocity-area methods. 

Unsteady flow models are also appropriate for evaluating the effects of changes in a managed flow regime on 
downstream conditions prior to the implementation of any changes. 

In the formulation of this standard considerable assistance has been derived from ISO/TR 1 1627 : 1998 
'Measurement of liquid flow in open channels — Computing stream flow using an unsteady flow model'. 

For the purpose of deciding whether a particular requirement of this standard is complied with, the final value, 
observed or calculated, expressing the result of a test or analysis shall be rounded off in accordance with 
IS 2 : 1960 'Rules for rounding off numerical values (revised)'. The number of significant places retained in 

the rounded off value should be the same as that of the specified value in this standard. 
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1 SCOPE 

This standard lays down a method for computing 
continuous records of stream flow in an open channel 
through the numerical solution of the one- 
dimensional unsteady flow equations. Such an 
approach is typically identified as an unsteady flow 
model and generally involves the use of computers 
for solution of the flow equations. 

The method described in this standard is applicable 
to steady and unsteady flows and to tidal flows in 
which there are no significant longitudinal and 
vertical density gradients. The method is considered 
equivalent to, or better than, the commonly used 
stage-fall-discharge technique [see 7.2 and Annex C 
of IS 15119 (Part 1)] because the method uses 
information on the physical characteristics of the 
channel, including the cross-sectional geometry, 
channel rugosity and channel slope, and the method 
is based on a mathematical description of the physics 
of fluid flow. 

This standard also describes the theoretical basis 
and fundamental assumptions of the technique, and 
provides a summary of selected numerical methods 
used to solve the unsteady flow equations. The 
details on the application of an unsteady flow model, 
including data requirements, procedures for model 
calibration, testing, and applications, and 
identification of uncertainties associated with the 
method have also been included. 

This standard, however does not provide sufficient 
information for the development of a computer 
programme for solving the unsteady flow equations, 
but rather it is based on the assumption that an 
adequately documented computer programme is 
available. 

2 REFERENCES 

The following standards contain provisions which 
through reference in this text, constitute provisions 
of this standard. At the time of publication, the 
editions indicated were valid. All standards are 
subject to revision and parties to agreements based 
on this standard are encouraged to investigate the 
possibility of applying the most recent editions of 
the standards indicated below: 



IS No. 



1191:2003 



1192:1981 



2912: 1999/ 
ISO 1070: 1992 

3910:1992 



3912:1993 

9163 

(Parti): 1979 



15118:2002/ 
IS04373: 1995 

15119 

(Parti): 2002/ 
ISO 1100-1:1996 

(Part2):2002/ 
ISO 1100-2: 1998 

15122:2002/ 
ISO 2425: 1999 

3 DEFINITIONS 



Tide 

Hydrometrie determinations — 
Vocabulary and symbols (second 
revision) 

Velocity area methods for 
measurement of flow in open 
channels {first revision) 

Liquid flow measurement in open 
channels — Slope — Area 
method {first revision) 

Requirements for rotating 
element current meters (cup type) 
for water flow measurement {first 
revision) 

Sounding rods — Functional 
requirements {first revision) 

Dilution methods of 

measurement of steady flow: 
Part 1 Constant rate injection 
method 

Measurement of liquid flow in 
open channel — Water level 
measuring devices 

Measurement of liquid flow in 
open channels: 

Establishment and operation of 
a gauging station 

Determination of the stage- 
discharge relation 

Measurement of liquid flow in 
open channels under tidal 
conditions 



For the purpose of this standard, the definitions 
given in IS 1 191 and the following shall apply. 

3.1 Boundary Condition — A boundary condition 
is a condition that a dependent variable of a 
differential equation must satisfy' along the boundary 
of the model domain. Boundary conditions for the 
dependent variables should be specified at the 
physical extremities of the modeled region for the 
duration of model application. 
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3.2 Courant Condition — The usual condition for 
the numerical stability of the explicit formulation of a 
numerical scheme which requires that the ratio of 
the propagation speed of a physical disturbance to 
that of a numerical signal should not exceed unity. 

3.3 Explicit Finite — Difference Numerical 
Scheme — Explicit numerical schemes convert either 
the characteristic equations or the governing 
equations to a system of linear algebraic equations 
from which the unknowns may be solved directly 
(explicitly) without iterative computations. 
Dependent variables on the advanced time level are 
determined one point at a time from known values 
and conditions at the present or previous time levels. 
Explicit schemes are only conditionally stable, 
meaning that errors may grow as the solution 
progresses, and the errors are a function of the time 
and distance finite-difference step sizes. Explicit 
schemes are generally stable when the courant 
conditions is met, which results in limitations on the 
distance step and maximum time which can be used. 

3.4 Gradually — Varied, Unsteady Flow — Generally 
non-uniform How in which there are no abrupt 
changes in depth along the longitudinal axis of the 
channel, and in which depth (and velocity and 
discharge) change with time. 

3.5 Hydrograph — A relation in graphical, 
equational, or tabular form between time and flow 
variables such as discharge depth, velocity, and 
stage. Stage and discharge hydrographs are typically 
used for open channel flows. 

3.6 Implicit Finite — Difference Numerical 

Scheme — Implicit numerical schemes convert either 
the characteristic equations or the governing 
equations to a system on non-linear algebraic 
equations from which the unknowns must be solved 
iteratively. All of the unknowns within the model 
domain are determined simultaneous!}, rather than 
point-by-point as with explicit methods. Implicit 
methods are generally stable, and are more 
computationally efficient than explicit schemes, but 
implicit schemes require more complex computer 
algorithms than do explicit schemes. 

3.7 Initial Conditions — A description of the 
dynamic conditions (typically, discharge and depth 
of flow for unsteady flow models) in the model 
domain at some specified time, usually the beginning 
of the simulation period. For ail subsequent times, 
the governing equations and the boundary- 
conditions describe the state of the system. 

3.8 Method of Characteristics — The method of 
characteristics is a mathematical approach for solving 
boundary-value problems by transforming the 
original partial differential equations representing 



the physical system into corresponding 
characteristic equations. The characteristic 
equations are ordinary differential equations and 
generally are more amenable to numerical solution 
than are the original partial differential equations. 

3.9 Momentum Coefficient — The momentum 
coefficient, also known as the Boussinesq coefficient, 
quantifies the deviation of the velocity at any point 
in a cross-section from uniform velocity distribution 
in the same cross-section. A value of unity indicates 
that a uniform velocity distribution is present in the 
cross-section. The momentum coefficient generally 
varies between about 1 .0 1 and 1 . 1 2 for fairly straight, 
prismatic channels; coefficients are typically smaller 
for large, deep channels than for small channels. 

4 UNITS OF MEASUREMENT 

The units of measurement used in this standard are 
SI units 

5 PRINCIPLES OF UNSTEADY FLOW MODELS 

5.1 Governing Equations 

The foundations for the fundamental derivation of 
the governing one-dimensional unsteady flow 
equations were laid by the 19th century hydraulicians 
Coriolis, Boussinesq, and Saint Venant. The 
governing equations are the one-dimensional, cross- 
sectionally averaged expressions for the following: 

a) Conservation of mass (or equation of 
continuity) 



8A dQ 

a & • 

b) Conservation of linear momentum 



..(1) 



ft M A 



where 



ax ' * - 



A = cross-sectional area of the channel, 

and varies with x, f, and r; 
t = time: 

O - discharge, and varies with x and t; 
u a longitudinal component of the 
lateral inflow velocity, and varies 
w ith x and /; 
x = longitudinal position along the 

channel axis; 
z = depth of flow, and varies w ith .v and /; 
g =■ acceleration of gravity; 
ft ** momentum coefficient, and varies 

with x T z, and /; 
q « lateral inflow per unit length of 
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channel, and varies with x and t; 
bed slope, and varies with x: and 
friction slope, and varies with x, / 
and z. The momentum coefficient 
may be computed as: 



pa" dA 



J U" A 
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where 



w ~ velocity in some elemental area 
dA 9 and 

U = mean velocity in the same 
cross-section having a total 
area A. 

The friction slope, 8 accounts for the resistance due 
to external boundary stresses. The friction slope is 
generally written as: 



«f = 



OlOtn' 



AR~ 



...(4) 



where 

R = hydraulic radius, and 
n = manning coefficient. 

Both R and n can vary as a function of x, z, and /. 
Equation 4 is based on the assumption that the 
manning equation for steady uniform flow provides 
a reasonable approximation for S f in unsteady, non- 
uniform flow. 

Equation 2 can be modified to include a term 

accounting for the momentum imparted to the water 
by a temporally and spatially varying wind. Equations 
1 and 2 also can be written with the following: 

a) Depth and velocity, 

b) Stage and velocity, or 

c) Stage and discharge as the dependent 
variables. 

Equations 1 and 2 apply to the unsteady, spatially- 
varied, turbulent free-surface flow of an 
incompressible viscous fluid in an open channel of 

arbitrary cross-section and alignment. The equations 
are solved simultaneously for the unknowns z (depth 
of flow) and Q (discharge) as a function of time (/) 

and longitudinal position (x). 

5.2 Assumptions Upon Which Governing Equations 
are Based 

Equations 1 and 2 are derived from first principles 
and may be obtained directly from the three- 
dimensional equation of mass continuity and the 
Navier-Stokes equations which are general, three- 
dimensional statements of the conservation of 
momentum for any fluid flow. A number of 
assumptions are required to derive Equations 1 and 



2. An unsteady flow model which is based on 
Equations 1 and 2 should generallv be applied to 
those conditions in which none of the major 
assumptions are severely violated. The assumptions 
are as follows: 

a) Flow is approximately one-dimensional, 
meaning that the predominant spatial variation 
in dynamic conditions (discharge, velocity and 
stage) is in the longitudinal direction. 

b) Fluid density is homogeneous throughout the 
modeled reach. 

c) Vertical accelerations are negligible (the 
hydrostatic pressure distribution is 
applicable). 

d) Velocity is uniformly distributed in a given 
cross-section. Inclusion of the momentum 
coefficient in Equation 2 allows this 
assumption to be violated somewhat, but there 
should be no flow separation, and streamline 
should not be highly curvilinear 

e) Neither aggradation nor degradation of the 
flow channel occurs. 

Turbulence and energy dissipation can be 
described by resistance laws formulated for 
steady, uniform flow (required for Equation 4). 

g) There are no abrupt changes in channel shape 
or alignment. The velocity is zero at the 
channel boundary. 

h) There is no superelevation of the water level 
at any cross-section. 

j) Surface tension and the density of air at the 
free surface are negligible. 

5.3 Simplified Models 

A number of techniques have been used to simplify 
Equations 1 and 2 to provide approximate unsteady 
flow models. These simplified models generally 
provide results with less computational effort and 
fewer data than is required for solution of the full 
equations. However, the models have limited 
applicability, and it is more appropriate to use a 
general unsteady flow model based on Equations 1 
and 2 to obtain reliable records of discharge under a 
wide range of conditions. A brief summary of 
simplified models follows. 

5.3.1 Empirical Models 

Empirical models are based on observations of past 
flood events. These models are limited to applications 
in which sufficient observations of inflows and 
outflows of a river section are available to calibrate 
essential empirical relations or routing coefficients. 
These models are typically applied to slowly 
fluctuating rivers with negligible lateral inflows and 
backwater effects. 
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5.3.2 Hydrotogic Models 

Hydrologic models are based on the continuity 
equation written as: 

I-O-dSdt -.-(5) 

where 

/ = inflow to the modelled river section, 
O - outflow from the section, and 
dS * change in storage within the section 
during the time interval dt 

The storage is generally assumed to be related to 
the inflow or outflow by some empirically-determined 
storage constant- Hvdrologic models are limited to 
applications, in which the stage-discharge relation 
is single-valued, and are not applicable to flows 
having backwater effects, significant lateral inflows, 
or looped stage-discharge relations. Difficulties in 
solving Equation 5 are often encountered when flows 
are changing rapidly with time. 

5.3.3 Linearized Models 

Linearized models are derived from Equations I and 
2 by ignoring or linearizing nonlinear terms in the 
equations. The linearized Equations can then be 
analytically integrated with less computational effort 
than is required for numerical integration of 
Equations 1 and 2. The most common simplifying 
assumptions for these models are: 

a) Acceleration term (second term) in the 
momentum Equation (equation 2) is 
negligible; 

b) Cross-sectional area A and channel bottom 
slope (S ) are constant; the friction slope (5 ) 
is linearized with respect to discharge and 
depth;. 

c) There is no lateral inflow; and 

d) Routed flood wave has a simple shape 
described by an analytical expression. 

These assumptions severely limit the applicability 
of linearized models. 

5.3.4 Kinematic Wave Model 

The kinematic wave model is derived by assuming 
that all terms in the momentum equation are negligible 
relative to the friction slope (S t ) and the bed slope 
(S ). and that there is no lateral inflow, so that 



S f S 
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As a consequence of Equation 6, the discharge for a 
kinematic flow is equal to the normal discharge. This 
means that the momentum of the unsteady flow is 
described by an expression such as the Manning or 
Chezy equations, in which flow is a single-valued 



function of depth of flow. Moreover, kinematic waves 
travel without attenuation of the peak flow, but the 
shape of the flood wave is modified as the wave is 
translated downstream. The kinematic wave model 
allows only the downstream propagation of flow 
disturbances so that backwater and tidal effects 
cannot be modeled. Numerous analytical solutions 
exist for applications of the kinematic wave model to 
specific flow geometries, and these models are most 
widely used in the routing of overland flow of 
precipitation runoff. 

5.3.5 Diffusion Analogy Model 

The diffusion analogy model is obtained by assuming 
that the channel is prismatic, that the local and 
convective acceleration terms in the momentum 
equation are negligible, and that there is no lateral 
inflow. The continuity and momentum equations may 
then be combined to form a single parabolic partial* 
differential equation, which is in the form of the so- 
called convective-diffusion equation with the single 
unknown of discharge. The local and convective 
acceleration terms, the first two terms in Equation 2, 
are often small in steep streams. 

The diffusion analogy model can be used to compute 
flows affected by backwater conditions. However, 
the diffusion model is limited to applications in which 
flows change relatively slowly, and in which the 
channel has a rather uniform geometry throughout 
the modeled reach. 

5.4 Numerical Techniques for Solution of 
Governing Equations 

No known analytical solutions exist for Equations 1 
and 2. Consequently, numerical techniques are used 
to convert Equations 1 and 2 into algebraic Equations 
that may be solved for z and Q at finite, incremental 
values ofx and /. This solution depends on the proper 
description of the cross-sectional area as a function 
of x and /, and on the availability of accurate 
boundary condition data. 

A variety of numerical techniques have been 
proposed and used to solve the unsteady flow 
equations. Although finite element methods may be 
used to solve the Equations, finite-difference 
techniques generally are more appropriate for the 
solution of the one-dimensional partial-differential 
Equations, describing unsteady open-channel flow. 
The three broad categories of numerical techniques 
are as given below: 

a) Method of characteristics, 

b) Explicit finite-difference methods, and 

c) Implicit finite-difference methods. 

Numerous variations of each of these general 
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categories of techniques exist. The methods are 
briefly reviewed to provide some perspective on 

advantages and disadvantages of each method. 

5.4. 1 Method of Characteristics 

The method of characteristics is a mathematical 
approach for solving boundary-value problems by 
transforming the original partial differential 
equations representing the physical system into 
corresponding characteristic equations. In this 
context, the characteristic is the speed of a wave 
relative to a stationary observer. Characteristic 
equations are ordinary differential equations and 
generally are more amenable to numerical solution 
than the original partial differential equations. The 
characteristic equations are solved using either 
explicit or implicit finite-difference methods. 

The method of characteristics can be used with a 
curvilinear grid or a rectangular grid in the x-t domain. 
The curvilinear grid general!) is not used for solution 
of the unsteady flow equations in natural open 
channels. The nature of characteristics is such that 
the wave trains in the x-t domain usually are not 
orthogonal, so solutions of the characteristic 
equations typically do not coincide with a point on 
the rectangular grid representing the natural system. 
Consequently, an interpolation scheme is required 
to transfer results from the characteristic network to 
the rectangular grid representing the flow system. 
The accuracy of the interpolation scheme plays a 
major role in determining the performance of the 
method of characteristics in solving the governing 
equations. 

5.4.2 Explicit Finite-difference Methods 

Explicit numerical schemes convert either the 
characteristic equations or the governing equations 
to a system of linear algebraic equations from which 
the unknowns may be solved directly (explicitly) 
without iterative computations. Dependent variables 
on the advanced time level are determined one point 
at a time from known values and conditions at the 
present or previous time levels. Explicit schemes are 
only conditionally stable, meaning that errors may 
grow as the solution progresses, and the errors are a 
function of the time and distance finite-difference 
step sizes. Explicit schemes are generally stable when 
the courant condition is met, which results in 
limitations on the distance step and maximum time 
which can be used. 

In order to meet numerical stability requirements, the 
computational time step must decrease as the 
hydraulic depth increases. Consequently, 
computational time steps may be required to be on 
the order of a few minutes for unsteady flow modeis 



of large rivers, which makes the models somewhat 
computationally inefficient. Explicit finite-difference 
schemes also require that the computational distance 
steps be equal throughout the model domain, which 
may be a disadvantage for some systems. 

5.4.3 Implicit Finite -difference Methods 

Implicit numerical schemes convert either the 

characteristic equations or the governing equations 
to a system on nonlinear algebraic equations from 
which the unknowns must be solved iteratively. 
Consequently, a system of 2 N algebraic equations 
is generated for a model having N cross-sections 
along the x-axis. All of the unknown within the model 
domain are determined simultaneously, rather than 
point-by-point as with explicit methods. 

Weighting factors are typically required in the 
application of implicit schemes. These factors 
determine the time between adjacent time levels at 
which (a) the spatial derivatives, and (b) functional 
quantities are evaluated; functional quantities are 
such features as cross-sectional area, top width, and 
hydraulic radius, all of which are functions of the 
computed depth of flow. Some judgement is required 
in selecting these weighting factors, and the 
weighting factors often are adjusted as part of the 
model calibration process. The accuracy of the 
numerical scheme generally decreases as the factor 
approaches one. where the terms in the governing 
equations are expressed entirely in terms of the future 
time step. 

Fewer numerical stability problems are encountered 
with implicit schemes than with explicit schemes. 
Numerical instabilities can occur when modeling 
rapidly varying flows if the time step is large and if 
the spatial derivatives are not sufficiently weighted 
toward the future time step. Non-linearities caused 
by irregular cross-sections having widths that vary 
rapidlv along the channel or with depth also can cause 
numerical instabilities in implicit models. 

6 DATA REQUIREMEiNTS 

Data are required to construct, calibrate, test, and 
applv unsteady flow models. Referenced Indian 
Standards for the measurement of velocity and 
discharge, for collection of water-level and discharge 
records should be followed. 

In general, data are required at model boundaries for 
the entire period for which flow is to be computed 
using the unsteady flow model. Short-term records 
and discrete measurements are needed at locations 
within the model domain for the period which is used 
for model calibration and testing. 
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6.1 Selection of Model Boundaries 

Reliable, accurate, and appropriate boundary 
condition data are required for successful 
computation of stream flow using an unsteady flow 
model. Model boundaries must be selected prior to 
the installation of data-collection instrumentation. 
Boundaries should be in locations where there are a 
minimum of flow disturbances, such as sharp bends, 
rapid changes in cross-sectional geometry, and major 
inflows. The modeled reach length also should be 
sufficiently long to permit accurate determination of 
the longitudinal water-surface slope so that adverse 
effects of measurement errors are minimized. 
Moreover, as subsequently discussed, discharge 
records should be used as the upstream boundary 
condition, if at all possible. For these reasons, it may 
be expedient to extend the model domain beyond the 
reach for which stream flow computations are actually 
needed in order to obtain the necessary data for 
model boundary conditions. 

6.2 Stage Data 

Stage data are required at all external boundaries of 
the modeled system in order to specify boundary 
conditions (Stage typically is required even if 
discharge is used as the upstream boundary 
condition because stage is generally needed in the 
computation of discharge whether from a stage- 
discharge relation or in-situ velocity meters). 
Multiple channel systems require stage data at 
upstream and downstream external boundaries of 
each channel. Stage data are not required at internal 
junctions where channels join within the model 
domain. Stage also should be measured in at least 
one, and preferably three, locations within the model 
domain to provide data for model calibration and 
testing. 

It is critically important that all stage measurements 
be referenced to a common datum. Errors in gauge 
datum translate into errors in water-surface , which 
greatly affects computed stream flow through the 
third term in Equation 2. Correct datums are 
particularly important (and perhaps more difficult to 
obtain) for low-gradient channels, where unsteady 
flow models are often applied. The use of discharge 
as the upstream boundary condition removes much 
of the uncertainty associated with potential errors in 
gauge datums. 

Except in very large rivers (widths of several hundred 
meters or more), stage should be measured at 1 5 min 
intervals or less for reliable modeling of flow 
transients. If possible, the stage measurement 
interval should be a whole multiple of the 
computational interval, and should be no more than 
about five times the computational interval. 



Synchronous measurement of stage at all recorders 
also is required for application of unsteady flow 
models for computing stream flow. Asynchronous 
measurements, like datum errors, translate into errors 
in water-surface slope and, hence, errors in computed 
stream flow. 

Stage should be measured following procedures 
outlined in IS 15119 (Part 1) using equipment 
described in IS 151 18 . 

6.3 Velocity Data 

Measurements of velocity in the study reach are 
required to (a) evaluate the assumption of one- 
dimensional flow, and (b) compute the momentum 
coefficient (Equation 3). Stream velocities obtained 
during discharge measurements generally are 
adequate for these purposes. 

Velocities should be measured following procedures 
outlined in IS 1 192 or IS 15122 using equipment 
described in IS 39 10 and IS 39 12. 

6.4 Discharge Data 

Discharge data are required for model calibration and 
may be needed as boundary data. Discharge data 
may be either (a) a continuous time series obtained 
from a stage-discharge rating or by using the in-situ 
velocity meters, such as ultrasonic velocity meters 
or (b) discrete measurements. Time-series of 
discharge are generally required only at the upstream 
boundaries. Discrete measurements are made within 
the model domain for the purposes of model 
calibration and testing. 

Discharge may be measured using methods 
described in IS 1192, IS 2912, IS 9163 (Part 1), 
IS 15119(Part 1), IS 151 19 (Part 2)and IS 15122. 

6.5 Lateral Inflows and Withdrawals 

Time-varying records of major inflows into and 
withdrawals from the modeled reach must be included 
in the model to maintain mass balance. Inflows or 
losses which are relatively constant throughout the 
modeled reach such as ground-water inputs or losses, 
can be lumped into a few discrete points in the reach 
or can be included at each computational node. 
Inflows from major streams must be gauged. Inflows 
from minor streams and local areas along the channel 
can be estimated using data from nearby gauged 
streams and drainage area ratios. 

6.6 Channel Cross-Section Data 

Cross-sectional data consists of a set of longitudinal, 
lateral and vertical coordinates which describe the 
location and configuration of the cross-section. A 
measured cross-section is typically required at each 



IS 15823: 2009 



computational node within the model. The exact 
spacing of the computational nodes, however, is not 
known until convergence testing has been completed 
and stability criteria have been evaluated. 

Numerical solutions of the governing equations 
generally use the average of the measured cross- 
sections at the upstream and downstream ends of a 
reach to represent the cross-sectional geometry of 
the entire reach bounded by the two measured 
sections. Consequently, measured cross-sections 
should be fairly representative of conditions 
upstream and downstream of the measurement. The 
measured sections also should be spaced sufficiently 
close so that large changes in channel geometry do 
not occur between the sections. 

The longitudinal distance between cross-sections 
can be determined from a map having a scale of 
1 :24 000 or larger. Longitudinal distances should be 
measured along the centerline of the channel. 

Cross-sections should be measured in accordance 
with IS 1 192 or by using a recording fathometer. It is 
generally not practical to run levels from an 
established benchmark to each cross-section. 
Consequently, cross-sections can be measured 
during steady flow when water levels are not 
changing, and the stage at the measurement site can 
be interpolated from measured stage at upstream and 
downstream stage recorders. Cross-sections shall be 
referenced to the same datum as the stage records. 



7 CALIBRATION AND TESTING 
UNSTEADY FLOW MODELS 



OF 



Modeling is based on abstraction of a physical 
system to a mathematical expression and replication 
of the system using these expressions and 
appropriate field data. The analyst should identify 
the important features of the flow system and ensure 
that those features are reflected in the model which 
is selected for application to the study reach. 
Important general model attributes include (a) the 
ability to simulate a wide range of flow conditions; 
(b) the ability to represent a range of complex channel 
conditions and geometric: (c) a stable, numerically- 
convergent, efficient computational scheme, and (d) 
a system for processing model input data and output 
simulation results. 

7.1 Preliminary Tests 

In many cases, it is appropriate to conduct 
preliminary tests using simplified channel geometry 
and boundary conditions with the unsteady flow 
model. Tests should be conducted if the model is 
poorly documented, or if the user is unfamiliar with 
the model. 

Tests should be conducted usine a channel which 



has a uniform, rectangular cross-section, and with 
the model configured for the study reach. Tests using 
the rectangular cross-section model could include 
the following: 

a) No inflow and no bed slope — no flow 
should be generated within the model 
domain. 

b) Steady inflow — mass should be conserved. 

c) Unsteady inflow — mass should be 
conserved. 

d) Triangular-shaped inflow hydrograph in 
channel with no bed slope-peak flow should 
not be significantly attenuated. 

After the model for the study reach has been 
constructed using the measured cross-sectional data, 
tests which should be performed include: 

a) No inflow or water surface slope — no flow 
should be generated in the model 
domain; this test also determines if 
unintentional openings in the boundaries 
are present. 

b) Steady and unsteady flow — mass should 
be conserved. 

c) Rapid change in inflow boundary 
conditions — no numerical instabilities 
should be generated. 

d) Change in boundary conditions from one 
steady flow to another flow — the amount 
of time required for all flows within the model 
domain to reach the new steady- 
state condition is an indication of how long 
initial conditions persist within the 
model domain; model results generally are 
not accepted until effects of the initial 
conditions are transported out of the model 
domain so that model is responding to 
boundary conditions only. 

Other tests may be performed as needed, but these 
simple tests can be used to document general model 
performance and should allow users to gain a better 
understanding of model capabilities and limitations 
prior to application to the study reach, 

7.2 Computational Grid and Time Step 

The computational grid is used to represent the 
physical system in the unsteady flow model. 
Junctions, inflows, and outflows should be 
represented by the computational grid. 

The channel system is subdivided into a number of 
finite segments for solution of the numerical 
approximation of the governing equations. The 
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solution points are either at the ends or the midpoint 
of each segment. The computational grid should be 
established such that computations of stage and 
discharge coincide with locations of data collection, 
or at locations where computed data are required. 
Measured cross-sectional data should be available 
at the ends of each computational segment or grid 
cell. 

Some models allow non-uniform segment lengths, 
but others require that all segments have the same 
length. Segment lengths should be at least three 
times greater than the width of the channel, and are 
frequently five to ten times the channel width. Exact 
segment length is determined during convergence 
testing. 

The representation of overbank areas, or the flood 
plain, in the model should be done with caution. For 
relatively narrow flood plains in which the flow length 
of the flood plain is very nearly equal to the flow 
length of the channel, model assumptions of one- 
dimensional flow are not violated. Broad flood plains, 
which bound a sinuous channel, may have distinctly 
different flow lengths than does the channel. 
Moreover, filling and draining of the floodpiain may 
lag the rise and fall of water in the channel. In these 
cases, the one-dimensional flow model described in 
this standard is not appropriate for application. 

In some instances, off-channel storage reservoirs 
have been included in one-dimensional flow models 
to account for storage and release of water from the 
flood plain. These formulations account for the slow 
storage and release of flood plain waters. However, 
such a model design also converts a process which 
occurs throughout some finite reach of the channel 
into a single point process. 

The computational time step should be sufficiently 
small to accurately represent flow transients which 
occur in the modeled system. Generally, the 
computational time step is reduced to meet stability 
criteria rather than adjusting the spatial discretization 
interval. As with the computational grid, convergence 
testing helps define the maximum time step which 
can be used. 

7.3 Convergence Testing 

A finite -difference solution to a partial-differential 
equation is spatially convergent, if the numerical 
solution approaches the true solution of the 
differential equation as the finite-difference spatial 
discretization approaches zero. Spatial convergent 
can be tested by repeatedly applying the model with 
a fixed set of boundary conditions for successively 
smaller computational discretizations. The model is 
spatially convergent if no further change in model 
results is observed as the spatial step is refined. 



Likewise, a model is temporally convergent, if model 
results remain substantially unchanged as the 
computational time step is decreased. To determine 
the effects of spatial discretization and time step on 
model results, convergence testing should be 
conducted prior to model calibration. 

7.4 Boundary and Initial Conditions 

Two initial conditions (Q and z for the formulation of 
the unsteady flow equations used in this standard) 
are required at each computational node in the model 
domain. For the initial application of the model to a 
study reach, common initial conditions are a steady 
flow, equal to the initial boundary-condition flow, 
and a water surface which slopes linearly from a 
measured upstream stage to a measured downstream 
stage. Output from a previous unsteady flow model 
application also may be used to determine initial 
conditions for a simulation which follows 
sequentially in time. 

The model may only be applied for periods which 
have measured boundary conditions. Boundary 
conditions include a time series of measured stage 
at the downstream boundaries, measured stage or 
discharge at the upstream boundaries, and measured 
lateral inflows. As previously discussed, better 
results often are obtained when discharge is used as 
the upstream boundary condition. 

7.5 Calibration 

Model calibration is required to adapt a general 
unsteady flow model to the specific application for 
which stream flows are to be computed. Calibrations 
are accomplished by adjusting model parameters 
until model results agree with observations. 

Essentially all components of the model are subject 
to adjustment during model calibration. Components 
that are directly measurable and physically well 
defined, however, are typically less subject to 
adjustment than are those that might not be directly 
measured. Measures for quantifying the calibration 
are given in 7.6. 

Initially, the resistance coefficient, momentum 
coefficient, and weighting coefficients for the 
numerical scheme should be varied because these 
parameters cannot be measured. Boundary gage 
datums may be adjusted slightly, if there is some 
uncertainty about the accuracy of the datum. Cross- 
sectional geometry also may be adjusted during the 
calibration process. Adjustment of channel geometry 
is justified because the measured cross-sections are 
used to represent the average conditions within a 
computational segment rather than the actual 
conditions at the measured cross-section. 
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It is entirely possible to achieve a well-calibrated 
model with empirical coefficients which bear little 
resemblance to those justified by the physics and 
setting of the study reach. In such a case, application 
of the model to another condition is not justified 
without another calibration. 

Because of the important assumptions of the 
governing equations is that the channel is stable, 
and most natural channels undergo continuous 
change, model calibrations should be repeated at 
least annually. The amount of data required for the 
subsequent model calibrations should not be as 
great as for the initial calibration, unless the user 
suspects dramatic changes in the stream geometry. 

7.6 Validation 

Validation refers to comparison of model results to 
measured stage and discharge. Calibration data 
should not be used for validation. 

Graphical comparisons of measured and simulated 

information are often used, but may be misleading. 
Quantitative measures of validation include measures 
of deviation and testing of the statistical significance 
of the deviation. Measures of deviation between 
model results and data include absolute and relative 
error, and root mean square deviations. 

Statistical tests of significance are necessary to 
determine whether deviations are meaningful or 
whether they are simply related to variability in data. 
However, statistically independent data are needed 
to test for significance. Consequently, data and model 
results must be sampled at intervals greater than the 
correlation time scale before applying certain 
statistical tests, such as the /-test. 

If possible, the model should be validated over the 
range of flow conditions for which the model is to be 
applied. The model may be applied to flows 
reasonably outside the range for which the model 
was tested, as long as conditions do not change 
appreciably. For example a model calibrated for flows 
less than bankfull should not be applied to simulate 
over bank flows. 

Greater emphasis should be placed in validating 
against discharge rather than stage, because the 
purpose of the model is to compute stream flow. 
However, in most case, stage data are much easier 
and less expensive to obtain than are discharge data, 
and more stage data will be available for use in 
calibration, validation, and testing. 

7.7 Sensitivity Testing 

Sensitivity testing consists of evaluating the 
sensitivity of model results to changes in selected 
model parameters. Parameters, or conditions, which 



should be included in the sensitivity testing include 
the resistance coefficient, the weighting coefficients 
used in the numerical scheme, the momentum 
coefficient, boundary gage datums (particularly, if 
stage is used as the upstream boundary condition), 
and channel geometry. The usual procedure is to 
successively increment the parameters by small 
amounts, apply the model, and compare the results 
with results from the calibrated model. A model which 
is highly sensitive to small changes in one or more 
parameters may become unstable for conditions 
outside those used for model calibration, and greater 
care must be taken when applying such a model. 

8 UNCERTAINTIES 

8.1 Governing Equations 

Equations 1 and 2 strictly apply to an inflnitesimally 
small volume at an instant in time. For the 
development and application of unsteady flow 
models, the equations are assumed to apply to some 
finite volume which may have a length on the order 
of hundreds of meters, a width of tens of metres, and 
a depth of several metres; the equations also are 
assumed to apply to some finite duration, which may 
be as much as 1 h . Uncertainties associated with the 
extrapolation of the differential equations to these 
finite volumes and times exist and are difficult to 
exactly quantify. However, the previously described 
convergence testing does provide an indication, 
along with the satisfaction of numerical stability 
criteria, that the appropriate time and distance 
discretization has been selected. 

Flow in open channels is three-dimensional in nature, 
but is approximated by the governing equations as 
having only variations in the longitudinal direction. 
In addition to assumption of one-dimensional flow, 
a number of other assumptions were used in the 
development of the unsteady flow equations. Some 
of these assumptions are restrictive and deviations 
from the assumption (such as no channel geometry 
change and no density gradients) will likely produce 
erroneous simulations. Other assumptions, such as 
negligible surface tension effects, will not be violated 
in open-channel flow. All deviations from model 
assumptions should be clearly documented, and 
model testing should be used to assist in quantifying 
the effects of the deviations on flow computations. 

8.2 Numerical Approximations to Governing 
Equations 

Equations 1 and 2 are simplified mathematical 
expressions for the complex three-dimensional 
turbulent flow field in an open channel. The 
expressions are further simplified by making 
numerical approximations of the partial differential 
equations. These approximations introduce further 
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uncertainty into the problem. The uncertainty, which 
is difficult to quantify for any given application, can 
be evaluated through testing and documentation of 
the numerical scheme. 

Clear documentation of the numerical scheme, 
including equations, discretization, and results of 
tests of the scheme, should be available. Certain 
numerical schemes, including variations of the 
method of characteristics, explicit, and implicit 
schemes for one-dimensional unsteady flow 
modeling, have general testing and documentation 
that is widely available. However, many schemes, 
which have been devised for specific problems, may 
require extensive evaluation by model users. 

Documentation of a numerical scheme should include 
a discussion of the relations among grid size, time 
step, and the stability and accuracy of the scheme. 
Numerical dispersion introduced by the scheme, 
which can result in the damping of a steep wave 
front, should be quantified. 

8 J Model Parameters 

A number of physical processes which are not 
explicitly expressed in the momentum equation 
(Equation 2) are combined into the channel rugosity 
coefficient (Manning n, in this case). The Manning 
n is an empirically derived coefficient which is an 
approximation of channel rugosity in a steady 
uniform flow. The parameter in Equation 2, however, 
is applied to unsteady flow, and, from the physical 
perspective, includes the effects of turbulent 
dissipation of energy by several processes. 

The Manning n value, as well as weighting factors 
used in certain numerical schemes cannot be 
measured, and these parameters are adjusted during 
model calibration to obtain agreement between model 



results and prototype measurements. Consequently, 
these model parameters also may include the effects 
of the following: 

a) Deviation of the modeled system from model 
assumptions, 

b) Effects of numerical approximations to the 
governing equations, and 

c) Errors in field measurements. 

The uncertainties associated with combining these 
processes into a few model parameters is difficult, if 
not impossible, to quantify. However, model 
parameters selected during the calibration and testing 
phase of model development should be within the 
range of previously published values. Otherwise the 
unrealistic model parameter vaJues may be masking 
some serious errors within the model-errors which 
may result in poor simulations of streamflow for 
conditions which were not evaluated during model 
testing. 

8.4 Data for Model Development, Testing and 
Application 

The collection of hydraulic field data is subject to a 
number of uncertainties. Quantification of 
uncertainties associated with the measurement of 
stage, velocity, discharge and cross-sectional data 
is discussed in previously published standards. 

Uncertainties associated with estimation of lateral 
inflows and withdrawals should be minimized by 
accurately measuring as many of the inputs and 
losses as possible. 

The effects of uncertainties associated with errors 
in gage datums can be minimized by using discharge 
as an upstream boundary condition rather than stage. 
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